1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.03
## LDevsum 1.01 1.03
## dh0 1.44 2.21
## dh1 5.72 10.41
## dh2 3.96 7.16
## dl0 1.31 1.75
## dl1 1.37 1.89
## dl2 1.25 1.63
## dl3 1.02 1.04
##
## Multivariate psrf
##
## 5.44
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1405.5374 | 24268.361 |
| DIC3 | 1243.8191 | 21376.472 |
| PWAIC | 89.9428 | 1524.853 |
| WAIC | 1312.6123 | 22101.590 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.02
## LDevsum 1.02 1.06
## dh0 4.37 7.72
## dh1 7.31 12.81
## dh2 3.59 6.71
## dl0 1.42 2.11
## dl1 1.42 2.07
## dl2 1.12 1.33
## dl3 2.24 3.69
##
## Multivariate psrf
##
## 7.27
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1377.23924 | 24118.439 |
| DIC3 | 1243.55301 | 21367.683 |
| PWAIC | 89.58375 | 1524.755 |
| WAIC | 1311.70933 | 22092.162 |